Army  Research  Laboratory 


Ab  Initio  Calculation  of 
Vibrational  Absorption  and 
Circular  Dichroism  Spectra  Using 
Density  Functional  Force  Fields 


P.  J.  Stephens 
F.  J.  Devlin 

UNIVERSITY  OF  SOUTHERN  CALIFORNIA 

C.  F.  Chabalowski 

U.S.  ARMY  RESEARCH  LABORATORY 

M.  J.  Frisch 

LORENTZIAN  INC. 


ARL-TR-695 


February  1995 


19950222  025 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  IS  UNLIMITED. 


NOTICES 


Destroy  this  report  when  it  is  no  longer  needed.  DO  NOT  return  it  to  the  originator. 


Additional  copies  of  this  report  may  be  obtained  from  the  National  Technical  Information 
Service,  U.S.  Department  of  Commerce,  5285  Port  Royal  Road,  Springfield,  VA  22161. 


The  findings  of  this  report  are  not  to  be  construed  as  an  official  Department  of  the  Army 
position,  unless  so  designated  by  other  authorized  documents. 


The  use  of  trade  names  or  manufacturers’  names  in  this  report  does  not  constitute 
endorsement  of  any  commercial  product. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  tor  this  collection  of  information  is  estimated  to  average  t  hour  per  response,  including  the  time  for  reviewing  instructions  searchingeristing  data  sources 
oatherinq  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information  Send  comments  regarding  this  burden  estimate  or 

collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services  Directorate  for  information  Operations  i and  Reports.  1215  Jefferson 
Davis  Highway.  Suite  1204.  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0188).  Washington.  DC  20503. 


1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

February  1995  _ Final,  January  1993-August  1994 


4.  TITLE  AND  SUBTITLE  I  5-  FUNDING  NUMBERS 


Ab  Initio  Calculation  of  Vibrational  Absorption  and  Circular  Dichroism  Spectra 
Using  Density  Functional  Force  Fields 


6.  AUTHOR(S) 

P.  J.  Stephens,*  F.  J.  Devlin,*  C.  F.  Chabalowski,  and  M.  J.  Frisch** 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Research  Laboratory 

ATTN:  AMSRL-WT-PC 

Aberdeen  Proving  Ground,  MD  21005-5066 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


PR:  1L161 102AH43 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


U.S.  Army  Research  Laboratory 

ATTN:  AMSRL-OP-AP-L 

Aberdeen  Proving  Ground,  MD  21005-5066 


11.  SUPPLEMENTARY  NOTES 

*P.  J.  Stevens  and  F.  J.  Develin  are  employed  by  the  University  of  Southern  California. 
”M.  J.  Frisch  is  employed  by  Lorentzian  Inc. 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

ARL-TR-695 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  is  unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

The  unpolarized  absorption  and  circular  dichroism  spectra  of  the  fundamental  vibrational  transitions  of  the  chiral 
molecule,  4-methyl-2-oxetanone,  are  calculated  ab  initio.  Harmonic  force  fields  are  obtained  using  Density  Functional 
Theory  (DFT),  MP2,  and  SCF  methodologies  and  a  [5s4p2d/3s2p]  (TZ2P)  basis  set.  DFT  calculations  use  the  Local  Spin 
Density  Approximation  (LSDA),  BLYP,  and  Becke3LYP  (B3LYP)  density  functionals.  Mid-IR  spectra  predicted  using 
LSDA,  BLYP,  and  B3LYP  force  fields  are  of  significantly  different  quality,  the  B3LYP  force  field  yielding  spectra  in 
clearly  superior,  and  overall  excellent,  agreement  with  experiment.  The  MP2  force  field  yields  spectra  in  slightly  worse 
agreement  with  experiment  than  the  B3LYP  force  field.  The  SCF  force  field  yields  spectra  in  poor  agreement  with 
experiment.  The  basis  set  dependence  of  B3LYP  force  fields  is  also  explored:  the  6-3 1G*  and  TZ2P  basis  sets  give  very 
similar  results  while  the  3-21G  basis  set  yields  spectra  in  substantially  worse  agreemepfr  with  experiment. 


14.  SUBJECT  TERMS 


vibration  circular  dichroism  spectra,  density  functional  theory,  MP2  force  fields 


15.  NUMBER  OF  PAGES 

24 


16.  PRICE  CODE 


17  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION  20.  LIMITATION  OF  ABSTRACT 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 


UNCLASSIFIED 


UNCLASSIFIED 


UNCLASSIFIED 


Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std  Z39-18 
298-102 


Intentionally  left  blank. 


< 


ii 


ACKNOWLEDGMENTS 


The  past  and  ongoing  support  of  studies  of  Vibrational  Circular  Dichroism  at  the  University  of 
Southern  California  (USC)  by  the  National  Science  Foundation  (NSF),  the  National  Institutes  of  Health 
(NIH),  the  North  Atlantic  Treaty  Organization  (NATO),  and  the  San  Diego  Supercomputer  Center  is 
gratefully  acknowledged.  Dr.  C.  F.  Chabalowski  would  like  to  thank  the  U.S.  Army  Edgewood  Research, 
Development,  and  Engineering  Center  (ERDEC)  for  partial  support  of  this  work. 


Aooegglon  Tor  , 


lins  SRAM 

MIC  SAB 

UnaaD.o«aee»d 

JgstlfLcatic 

a 

□ 

Availability 

list 

V1 

Ivail  i 
Spec, 

and/cs* 

tal 

Intentionally  left  blank. 


iv 


TABLE  OF  CONTENTS 


Page 


ACKNOWLEDGMENTS  .  iii 

LIST  OF  FIGURES  .  vii 

1.  INTRODUCTION  . 1 

2.  METHODS .  2 

3.  RESULTS  AND  DISCUSSION .  4 

4.  CONCLUSION  .  9 

5.  REFERENCES .  13 

DISTRIBUTION  LIST .  15 


V 


Intentionally  left  blank. 


LIST  OF  FIGURES 


(a)  Experimental  unpolarized  absorption  spectrum  of  1  (in  CG4  and  CS2  solutions, 
from  Devlin  et  al.  [to  be  published]); 

(b)  Fundamental  unpolarized  absorption  spectrum  of  1  (obtained  from  (a)  by 
Lorentzian  fitting  [Devlin  et  al.,  to  be  published]); 

(cXg)  Calculated  unpolarized  absorption  spectrum  of  1  (In  (c),  (d),  and  (e),  force 
fields  and  APTs  are  calculated  using  DFT  and  B3LYP,  BLYP,  and  LSDA 
density  functionals,  respectively.  In  (f)  and  (g),  force  fields  and  APTs 
are  calculated  using  MP2  and  SCF  methodologies.  All  calculations  use 
the  TZ2P  basis  set.  Lorentzian  band  half-widths,  y,  are  4.0  cm'1.) . 

(a)  Experimental  VCD  spectrum  of  R-l  (in  CC14  and  CS2  solutions,  from 
Devlin  et  al.  [to  be  published]); 

(b)  Fundamental  VCD  spectrum  of  R-l  (obtained  from  (a)  by  Lorentzian  fitting 
[Devlin  et  al.,  to  be  published]); 

(cHg)  Calculated  VCD  spectrum  of  R-l  (These  originate  as  in  Figure  1.  All 
calculations  use  SCF  TZ2P  local  AATs.  Lorentzian  band  half-widths, 
y,  are  4.0  cm-1.)  . 

(a)  Fundamental  unpolarized  absorption  spectrum  of  1  (as  in  Figure  lb); 

(b) -(d)  Calculated  unpolarized  absorption  spectrum  of  1  (In  (b),  (c),  and  (d), 

force  fields  and  APTs  are  calculated  using  DFT,  the  B3LYP  density 
functional  and,  respectively,  the  TZ2P,  6-3 1G*  and  3-21G  Basis  sets. 
Lorentzian  band  half-widths,  y,  are  4.0  cm-1.) . . 

(a)  Fundamental  VCD  spectrum  of  R-l  (as  in  Figure  2b); 

(b) -(d)  Calculated  VCD  spectrum  of  R-l  (These  originate  as  in  Figure  3.  All 

calculations  use  SCF  TZ2P  local  AATs.  Lorentzian  band  half-widths, 
y,  are  4.0  cm-1.)  . 


Intentionally  left  blank. 


1.  INTRODUCTION 


The  use  of  density  functional  theory  (DFT)  in  the  ab  initio  calculation  of  molecular  properties  has 
recently  increased  dramatically  (Ziegler  1991;  Labanowski  and  Andzelm  [editors]  1991).  This  can  be 
attributed  to  the  following:  (1)  the  development  of  new  and  more  accurate  density  functionals;  (2)  the 
increasing  versatility,  efficiency,  and  availability  of  DFT  codes;  and,  most  importantly,  (3)  the  superior 
ratio  of  accuracy  to  effort  exhibited  by  DFT  computations  relative  to  other  ab  initio  methodologies. 

We  focus  here  on  the  application  of  DFT  to  the  prediction  of  vibrational  spectra.  With  the  most 
accurate  density  functionals  currently  in  use,  DFT  predicts  harmonic  vibrational  frequencies  of 
substantially  higher  accuracy  than  obtained  via  SCF  calculations  and  of  similar  accuracy  to  the  predictions 
of  MP2  calculations  (Johnson,  Gill,  and  Pople  1993;  Handy,  Murray,  and  Amos  1993).  Since  DFT 
calculations  are  substantially  less  demanding  computationally  than  MP2  calculations,  it  is  clear  that  DFT 
offers  significant  advantages  in  predicting  harmonic  vibrational  force  fields,  frequencies,  and  spectra. 

In  particular,  DFT  promises  to  greatly  facilitate  the  calculation  of  Vibrational  Grcular  Dichroism 
(VCD)  spectra  (Stephens  and  Lowe  1985).  Experience  has  shown  that  (with  the  exception  of  small, 
symmetrical  molecules)  the  SCF  harmonic  force  fields  are  generally  of  insufficient  accuracy  to  provide 
useful  predictions  of  VCD  spectra  (Kawiecki  1988;  Jalkanen  1989;  Bursi  1991).  Calculations  using  MP2 
harmonic  force  fields  are  substantially  more  successful  (Bursi  1991;  Amos,  Handy,  and  Palmieri  1990; 
Stephens  et  al.  1993,  Stephens  et  al.,  to  be  published;  Devlin  and  Stephens  1994),  but  also  much  more 
demanding.  The  lower  computational  demands  of  DFT  suggest  that  DFT  may  become  preferred  over  the 
MP2  methodology  in  the  future  ab  initio  predictions  of  VCD  spectra. 

We  report  here  calculations  of  the  unpolarized  absorption  and  circular  dichroism  spectra  arising  from 
the  fundamental  vibrational  excitations  of  the  chiral  molecule  4-methyl-2-oxetanone,  1,  previously  studied 
(both  experimentally  and  theoretically)  in  our  laboratories  (Jalkanen  1989;  Devlin  et  al.,  to  be  published). 
We  compare  the  accuracies  of  spectra  predicted  using  DFT,  MP2,  and  SCF  harmonic  force  fields.  The 
accuracy  of  DFT  calculations  is  dependent  on  the  density  functional  employed.  In  this  work,  we  compare 
three  significantly  different  density  functionals:  Local  Spin  Density  Approximation  (LSDA),  BLYP,  and 
Becke3LYP.  The  LSDA  and  BLYP  functionals  have  been  widely  used  (Ziegler  1991;  Labanowski  and 
Andzelm  [editors]  1991;  Johnson,  Gill,  and  Pople  1993;  Handy,  Murray,  and  Amos  1993).  The 
Becke3LYP  functional  is  a  hybrid  of  several  components,  whose  relative  weights  are  chosen  by  reference 
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to  experimental  thermochemical  data.  The  accuracies  of  such  hybrid  functionals  in  predicting  molecular 
geometries  and  vibrational  frequencies  have  not  yet  been  thoroughly  characterized. 

Our  calculations  take  advantage  of  recent  developments  in  the  computation  of  DFT  and  MP2  harmonic 
force  fields  via  analytical  derivative  techniques  (Trucks  et  al.,  unpublished;  Johnson  and  Frisch  1993, 
1994). 

2.  METHODS 


Harmonic  force  fields  for  1  were  calculated  ab  initio  at  SCF  and  MP2  levels  and  using  DFT  via 
GAUSSIAN  92  and  GAUSSIAN  92/DFT  (Frisch  et  al.  1992, 1993).  All  force  fields  were  calculated  using 
analytical  derivative  methods  (Trucks  et  al.,  unpublished;  Johnson  and  Frisch  1993,  1994).  DFT 
calculations  were  carried  out  using  three  density  functionals: 


(1)  LSDA  (Local  Spin  Density  Approximation).  This  uses  the  standard  local  exchange  functional 
(Becke  1989)  and  the  local  correlation  functional  of  Vosko,  Wilk,  and  Nusair  (VWN)  (1980). 

(2)  BLYP.  This  combines  the  standard  local  exchange  functional  with  the  gradient  correction  of 
Becke  (1989)  and  uses  the  Lee-Yang-Parr  (1988)  correlation  functional  (which  also  includes 
density  gradient  terms). 

(3)  Becke3LYP.  This  functional  is  a  hybrid  of  exact  (Hartree-Fock)  exchange  with  local  and 
gradient-corrected  exchange  and  correlation  terms,  as  first  suggested  by  Becke  (1993).  The 
exchange-correlation  functional  proposed  and  tested  by  Becke  was: 


Exc  -  0  - 


,  —  LSDA 
ao)Ev  +  a 


r-HF  .  _B88  —LSDA 
0EX  +  axAEx  +  Ec 


+  acAEc 


PW91 


(1) 


Here  AEX  8  is  Becke’s  gradient  correction  to  the  exchange  functional  and  AE^W91  is  the  Perdew-Wang 
gradient  correction  to  the  correlation  functional  (Perdew  1991).  Becke  suggested  coefficients  ao  =  0.2, 
ax  =  0.72,  and  ac  =  0.81,  based  on  fitting  to  heats  of  formation  of  small  molecules.  Only  single-point 
energies  were  involved  in  the  fit;  no  molecular  geometries  or  frequencies  were  used.  The  Becke3LYP 
functional  in  Gaussian  92/DFT  uses  the  values  of  a^  ax,  and  ac  suggested  by  Becke,  but  uses  LYP  for  the 
correlation  functional.  Since  LYP  does  not  have  an  easily  separable  local  component,  the  VWN  local 
correlation  expression  has  been  used  to  provide  the  different  coefficients  of  local  and  gradient  corrected 
correlation  functionals: 
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,  B3LYP 
'xc 


, —LSDA  _HF  .  CB88  „  _LYP  /t  „  ,CVWN  /0s 

(1  -  ao)Ex  +  aoEx  +  axAEx  +  acEc  +  (1  -  ac)Ec  .  (2) 


The  standard  fine  grid  in  Gaussian  92/DFT  (Trucks  and  Frisch,  to  be  published)  was  used  in  all  DFT 
calculations.  This  grid  was  produced  from  a  basic  grid  having  75  radial  shells  and  302  angular  points  per 
radial  shell,  leaving  about  7,000  points  per  atom  while  retaining  similar  accuracy  to  the  original  (75,302) 
grid.  Becke’s  (1988)  numerical  integration  techniques  were  employed. 

Atomic  Polar  Tensors  (APTs)  (Stephens  et  al.  1990)  were  calculated  at  the  same  time  as  harmonic 
force  fields.  Atomic  Axial  Tensors  (AATs)  were  calculated  using  the  Distributed  Origin  gauge  (Stephens 
et  al.  1990;  Stephens  1987),  in  which  the  AAT  of  nucleus  X,  (M^)0,  with  respect  to  origin  O  is  given 
by 


0< 


)°- 


—  i 

4hc 


ePr8  RXy  Pa5 


(3) 


where  R®  is  the  equilibrium  position  of  nucleus  X  relative  to  origin  O,  P^p  is  the  APT  of  nucleus  X,  and  (I^p)*' 
is  the  electronic  AAT  of  nucleus  X  calculated  with  the  origin  at  .  "Local"  AATs,  (Iap)\  are  calculated 
using  CADPAC  (version  5)  (Amos  1993)  at  the  SCF  level.  (At  the  present  time,  neither  MP2  nor  DFT 
codes  for  AATs  are  available). 

Vibrational  frequencies,  dipole  strengths,  and  rotational  strengths  were  calculated  from  the  harmonic 
force  fields,  APTs,  and  AATs  (Stephens  1985;  Jalkanen  et  al.  1990;  Kawiecki  et  al.  1991).  Unpolarized 
absorption  and  circular  dichroism  spectra  were  synthesized  thence  using  Lorentzian  band  shapes  (Kawiecki 
et  al.  1988).  Calculations  used  3-21G,  6-31G*  (Hehre  et  al.  1986)  and  [5s4p2d/3s2p],  TZ2P  (Stephens 
et  al.  1990)  basis  sets  (66,  102,  and  228  basis  functions,  respectively).  SCF  and  MP2  calculations  were 
earned  out  using  the  San  Diego  Supercomputer  Center  CRAY-YMP  and  C90  machines.  DFT  calculations 
(  were  carried  out  at  Lorentzian  Inc.  using  IBM  RS/6000-590  and  Silicon  Graphics  Incorporated  (SGI) 

Challenge  machines. 
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3.  RESULTS  AND  DISCUSSION 


Unpolarized  absorption  and  circular  dichroisra  spectra  of  1,  measured  in  CS2  and  CC14  solutions 
(Devlin  et  al.,  to  be  published),  are  reproduced  in  Figures  1  and  2  over  the  spectral  range  650-1,500  cm-1. 
Frequencies,  dipole  strengths,  and  rotational  strengths  of  transitions  assigned  to  fundamental  excitations, 
obtained  by  Lorentzian  fitting  of  the  experimental  spectra  (Devlin  et  al.,  to  be  published),  are  given  in 
Table  1.  The  spectra  attributable  to  these  transitions  alone  are  also  shown  in  Figures  1  and  2. 


Table  1.  Frequencies,  Dipole  Strengths,  and  Rotational  Strengths  of  la 


MP2 

j  DFT/B3LYP 

I  EXPl^ 

V 

D 

R 

V 

D 

R 

V 

D 

R 

1,528 

14 

-5 

1,500 

13 

-5 

1,453 

25 

-2 

1,516 

8 

4 

1,489 

8 

5 

1,441 

3 

3 

1,479 

47 

-11 

1,462 

43 

-11 

1,419 

35 

-9 

1,438 

52 

11 

1,422 

51 

11 

1,387 

58 

13 

1,394 

52 

27 

1,389 

50 

28 

1,350 

40 

42 

1,320 

97 

9 

1,316 

112 

3 

1,284 

161 

26 

1,230 

47 

-11 

1,229 

14 

-4 

1,198 

47 

-21 

1,216 

14 

0 

1,204 

19 

-5 

1,178 

35 

-17 

1,146 

191 

82 

1,137 

327 

99 

1,118 

483 

178 

1,125 

219 

-40 

1,113 

87 

-58 

1,099 

51 

-83 

1,077 

10 

-12 

1,076 

11 

-6 

1,055 

40 

-23 

1,056 

271 

-50 

1,041 

268 

-50 

1,022 

251 

-102 

973 

126 

10 

966 

91 

7 

959 

70 

12 

914 

9 

-29 

914 

8 

-27 

896 

10 

-54 

873 

221 

33 

850 

310 

46 

836 

371 

78 

829 

97 

26 

817 

34 

10 

812 

55 

29 

714 

3 

-3 

712 

2 

2 

711 

13 

15 

a  Frequencies,  v,  in  cm-1;  dipole  strengths  D  in  10'4Oesu2cm2;  rotational  strengths  R  in  KT^es^cm2 
R  values  are  for  R-(+)-l. 

Devlin  et  alM  to  be  published. 


4 


6 


Harmonic  frequencies,  dipole  strengths,  and  rotational  strengths  have  been  calculated  for  1  using  DFT 
at  the  TZ2P  basis  set  level  and  utilizing  the  LSDA,  BLYP,  and  B3LYP  density  functionals.  Unpolarized 
absorption  spectra  and  VCD  spectra  derived  thence  are  displayed  in  Figures  1  and  2.  The  absorption  and 
VCD  spectra  predicted  using  the  LSDA,  BLYP,  and  B3LYP  functionals  are  substantially  different, 
demonstrating  that  the  choice  of  functional  is  of  importance.  Comparison  to  the  fundamental  absorption 
and  VCD  spectra  shows  that  the  B3LYP  functional  provides  spectra  in  best  agreement  with  experiment 
when  judged  on  the  basis  of  the  qualitative  pattern  of  frequencies  and  intensities.  The  BLYP  and  LSDA 
functionals  yield  spectra  in  significantly  worse  agreement.  The  variations  among  predicted  spectra 
originate  predominantly  in  the  variations  in  the  harmonic  force  fields  (variations  in  the  APTs  are  of  minor 
consequence  in  comparison).  We  therefore  conclude  that  the  B3LYP,  BLYP,  and  LSDA  force  fields  are 
of  significantly  different  accuracy,  the  B3LYP  force  field  being  the  most  accurate. 

The  qualitative  agreement  of  the  absorption  and  VCD  spectra  calculated  using  the  B3LYP  functional 
with  the  experimental  spectra  is  impressive.  The  quantitative  agreement  of  theory  and  experiment  can  be 
assessed  from  the  calculated  and  experimental  frequencies,  dipole  strengths,  and  rotational  throughout,  the 
average  and  maximum  deviations  being  25  cm-1/2.0%  and  48  cm_1/3.3%.  Some  fraction  of  the 
differences  in  calculated  and  experimental  frequencies  can,  of  course,  be  attributed  to  anhaimonicity  and 
solvent  effects.  The  calculated  dipole  strengths  are  in  excellent  qualitative  agreement  with  experimental 
values;  quantitative  differences  are  attributable  to  a  combination  of  residual  errors  in  the 
DFT/B3LYP/TZ2P  calculation,  anharmonicity  and  solvent  effects,  and  to  experimental  error.  Calculated 
rotational  strengths  are  also  in  excellent  qualitative  agreement  with  experimental  values.  The  quantitative 
agreement  is  somewhat  poorer  than  in  the  case  of  dipole  strengths.  This  is  undoubtedly  attributable  to 
the  lower  accuracy  of  the  AATs,  compared  to  the  APTs,  due  to  the  absence  of  correlation  in  the 
calculation  of  local  AATs. 

The  results  of  TZ2P  MP2  calculations  (Chabalowski,  Devlin,  and  Stephens  to  be  submitted)  are  also 
given  in  Table  1  and  in  Figures  1  and  2.  The  agreement  with  experiment  of  the  spectra  calculated  using 
the  MP2  force  field  is  slightly  worse  than  obtained  from  the  B3LYP  force  field.  The  most  prominent 
difference  is  in  the  relative  absorption  intensities  of  the  bands  calculated  between  1,100  and  1,150  cm-1. 
On  the  other  hand,  the  agreement  is  significantly  better  than  obtained  from  the  BLYP  and  LSDA  force 
fields.  We  conclude  that  the  MP2  force  field  is  respectively  lower,  higher,  and  higher  in  accuracy  than 
the  B3LYP,  BLYP,  and  LSDA  force  fields. 
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The  quantitative  agreement  of  frequencies,  dipole  strengths,  and  rotational  strengths  calculated  from 
the  MP2  force  field  with  experimental  values  can  be  assessed  from  Table  1,  Calculated  frequencies  are 
in  all  cases  greater  than  experimental  frequencies,  the  average  and  maximum  deviations  being 
36  cm-1/3.0%  and  75  cm“V5.2%.  The  differences  are  in  all  cases  greater  than  for  the  B3LYP  force  field, 
and,  consequently,  so  are  the  average  and  maximum  deviations.  Calculated  dipole  strengths  are,  in  most 
cases,  further  from  experimental  values  than  those  obtained  from  the  B3LYP  force  field.  Comparison  of 
both  frequencies  and  dipole  strengths  thus  support  the  conclusion  that  the  harmonic  force  field  is  less 
accurately  calculated  by  the  MP2  methodology.  Rotational  strengths  calculated  from  the  MP2  force  field 
are  in  good  qualitative  agreement  with  both  experimental  values  and  the  predictions  of  the  B3LYP  force 
field.  Neither  of  the  calculations  exhibits  clear-cut  quantitative  superiority. 

Absorption  and  VCD  spectra  obtained  from  the  SCF  harmonic  force  field  are  also  shown  in  Figures  1 
and  2.  Both  absorption  and  VCD  spectra  are  qualitatively  different  from  the  experimental  spectra,  the 
differences  being  sufficiently  large  that  a  one-to-one  correspondence  between  calculated  and  observed 
bands  cannot  be  arrived  at;  i.e.,  the  predictions  are  insufficiently  accurate  to  yield  a  convincing  assignment 
of  the  experimental  spectra.  It  is  obvious  that  the  SCF  force  field  is  drastically  lower  in  accuracy  than 
the  DFT  and  MP2  force  fields. 

The  calculations  presented  in  Figures  1  and  2  and  in  Table  1  use  a  large  basis  set — TZ2P — in  order 
to  reduce  basis  set  errors  to  the  maximum  extent.  Having  established  the  superiority  of  the  B3LYP 
functional  over  the  LSDA  and  BLYP  functionals,  it  is  of  interest  to  examine  the  basis  set  dependence  of 
the  spectra  predicted  using  the  B3LYP  functional.  Figures  3  and  4  compare  absorption  and  VCD  spectra 
calculated  using  the  3-2 1G  and  6-3 1G*  basis  sets  to  those  obtained  at  the  TZ2P  level  and  to  experiment. 
Local  AATs  are  calculated  at  the  TZ2P  basis  set  level  throughout,  and  the  variations  in  spectra  are  due 
entirely  to  variations  in  the  DFT  force  fields  and  APTs  with  basis  sets.  The  absorption  and  VCD  spectra 
predicted  using  the  6-3 1G*  and  TZ2P  basis  sets  are  identical  qualitatively  and  very  little  different 
quantitatively.  The  agreement  with  the  experimental  spectra  is  almost  equally  good.  The  3-21G  basis  set 
yields  spectra  substantially  different  from  those  obtained  using  the  6-3 1G*  and  TZ2P  basis  sets  and  in 
much  worse  agreement  with  experiment. 
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4.  CONCLUSION 


In  the  case  of  1,  and  the  TZ2P  basis  set  level,  the  DFT  force  field  calculated  using  the  B3LYP  density 
functional  yields  mid-infrared  (IR)  absorption  and  VCD  spectra  in  impressive  agreement  with  experiment. 
DFT  force  fields  calculated  using  the  LSD  A  and  BLYP  density  functionals  yield  spectra  of  significantly 
lower  accuracy.  The  MP2  and  SCF  force  fields  yield  spectra  slightly  less  and  much  less  accurate, 
respectively,  than  the  DFT/B3LYP  force  field.  Calculations  on  a  wide  variety  of  molecules  are  clearly 
required  to  define  the  generality  of  these  results.  Such  calculations  are  under  way. 

While  the  LSDA  and  BLYP  functionals  have  been  widely  used  in  DFT  calculations  of  vibrational 
frequencies  (Johnson,  Gill,  and  Pople  1993;  Handy,  Murray,  and  Amos  1993),  this  is  not  the  case  for  the 
hybrid  functionals  recently  introduced  by  Becke  (1993).  Since  these  functionals  incorporate  parameters 
determined  by  fitting  to  experimental  thermochemical  data,  their  accuracy  in  predicting  molecular 
geometries  and  vibrational  force  fields  was  not  obvious.  Our  results  for  1  using  the  B3LYP  functional 
support  the  conclusion  that  the  methodology  of  Becke  leads  to  molecular  force  fields,  as  well  as 
thermochemical  properties,  of  improved  accuracy.  Our  results  also  show  that  B3LYP  calculations 
converge  rapidly  with  increasing  basis  set  size  and  that  the  cost-to-benefit  ratio  is  optimal  at  the  6-3 1G* 
basis  set  level.  6-3 1G*  will  be  the  basis  set  of  choice  in  B3LYP  calculations  on  much  larger  molecules. 

Our  calculations  clearly  demonstrate  the  sensitivity  of  predictions  of  vibrational  spectra  to  the  accuracy 
of  the  force  field,  and,  in  the  case  of  DFT  calculations,  to  the  choice  of  density  functional.  The 
comparison  to  experimental  spectra  of  vibrational  spectra  predicted  from  future  generations  of  density 
functionals  should  be  of  significant  utility  in  defining  their  relative  accuracies.  In  the  case  of  chiral 
molecules,  the  information  content  of  vibrational  spectroscopy  is  considerably  enhanced  when  absorption 
and  VCD  spectra  are  utilized  in  combination. 

VCD  spectroscopy  is  potentially  a  powerful  tool  in  elucidating  the  stereochemistry  of  chiral  molecules. 
The  work  reported  here  clearly  supports  the  expectation  that  DFT  now  provides  a  basis  for  calculations 
of  VCD  spectra  at  a  useful  level  of  accuracy  for  molecules  much  larger  than  1.  This  capability  will  be 
documented  in  forthcoming  publications. 
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wavenumbers 


Figure  3.  (a)  Fundamental  unpolarized  absorption  spectrum  of  1  (as  in  Figure  lb): 

(b)-(d)  Calculated  unpolarized  absorption  spectrum  of  1  (In  (b),  (cl,  and  (d).  force  fields  and 
APTs  are  calculated  using  DFT,  the  B3LYP  density  functional,  and,  respectively,  the 
TZ2P,  6-31G*.  and  3-21G  Basis  sets.  Lorentzian  band  half-widths,  v.  are  4.0  cm"1.). 
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Figure  4.  (a)  Fundamental  VCD  spectrum  of  R-l  (as  in  Figure  2b); 

(b)-(d)  Calculated  VCD  spectrum  of  R-l  (These  originate  as  in  Figure  3.  All  calculations  use 
SCF  TZ2P  local  AATs.  Lorentzian  band  half-widths,  y,  are  4.0  cm-1.). 
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DEPT  OF  CHEMISTRY 
BAKER  LAB 
ITHACA  NY  14853 

1  UNIV  OF  DELAWARE 

ATTN  T  BREL 
CHEMISTRY  DEPT 
NEWARK  DE  19711 

1  UNIV  OF  FLORIDA 

ATTN  J  WINEFORDNER 
DEPT  OF  CHEMISTRY 
GAINESVILLE  FL  32611 

3  GEORGIA  INSTITUTE  OF  TECHNOLOGY 

ATTN  E  PRICE 
W  C  STRAHLE 
BTZINN 

SCHOOL  OF  AEROSPACE  ENGRG 
ATLANTA  GA  30332 

1  UNIV  OF  ELINOIS 

ATTN  H  KREER 
DEPT  OF  MECH  ENGRG 
144MEB  1206  W  GREEN  ST 
URBANA  E  61801 


1  THE  JOHNS  HOPKINS  UNIV 

ATTN  T  W  CHRISTIAN 
CHEMICAL  PROPULSION 
INFORMATION  AGENCY 
10630  LITTLE  PATUXENT  PKWY 
STE  202 

COLUMBIA  MD  21044  3200 

1  UNIV  OF  MICHIGAN 

ATTN  G  M  FAETH 
GAS  DYNAMICS  LAB 
AEROSPACE  ENGRG  BLDG 
ANN  ARBOR  MI  48109  2140 

1  UNIV  OF  MINNESOTA 

ATTN  E  FLETCHER 
DEPT  OF  MECHANICAL  ENGRG 
MINNEAPOLIS  MN  55455 

4  PENNSYLVANIA  STATE  UNIV 

ATTN  K  KUO 
MMICCI 
S  THYNELL 
V  YANG 

DEPT  OF  MECHANICAL  ENGRG 
UNIV  PK  PA  16802 

1  POLYTECHNIC  INSTITUTE  OF  NY 
ATTN  S  LEDERMAN 
GRADUATE  CTR 

RTE  110 

FARMINGDALE  NY  11735 

2  PRINCETON  UNIV 
ATTN  K  BREZINSKY 
I  GLASSMAN 
FORRESTAL  CAMPUS  LIB 
PO  BOX  710 
PRINCETON  NJ  08540 

1  PURDUE  UNIV 

ATTN  J  R  OSBORN 

SCHOOL  OF  AERONAUTICS  &  ASTRONAUTICS 

GRISSOM  HALL 

WEST  LAFAYETTE  IN  47906 

1  PURDUE  UNIV 

ATTN  E  GRANT 
DEPT  OF  CHEMISTRY 
WEST  LAFAYETTE  IN  47906 
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2  PURDUE  UNIV 

ATTN  N  M  LAURENDEAU 
SNBMURTHY 

SCHOOL  OF  MECHANICAL  ENGRG 
TSPC  CHAFFEE  HALL 
WEST  LAFAYETTE  IN  47906 

1  RENSSELAER  POLYTECHNIC  INST 
ATTN  AFONTDN 
DEPT  OF  CHEMICAL  ENGRG 
TROY  NY  12181 

1  STANFORD  UNIV 

ATTN  R  HANSON 
DEPT  OF  MECHANICAL  ENGRG 
STANFORD  CA  94305 

1  UNIV  OF  TEXAS 

ATTN  W  GARDINER 
DEPT  OF  CHEMISTRY 
AUSTIN  TX  78712 

1  VA  POLYTECHNIC  INST  &  STATE  UNIV 

ATTN  J  A  SCHETZ 
BLACKSBURG  VA  24061 

1  APPLIED  COMBUSTION  TECH  INC 

ATTN  A  M  VARNEY 

PO  BOX  607885 
ORLANDO  FL  32860 

2  APPLIED  MECHANICS  REVIEWS 
ATTN  R  E  WHITE 

A  B  WENZEL 

THE  AMERICAN  SOCIETY  OF  MECH  ENGRG 

345  E  47TH  ST 

NEW  YORK  NY  10017 

1  ATLANTIC  RSCH  CORP 

ATTN  R  H  W  WAESCHE 
7511  WELLINGTON  RD 
GAINESVILLE  VA  22065 

1  TEXTRON  DEFENSE  SYSTEMS 

ATTN  A  PATRICK 
2385  REVERE  BEACH  PKWY 
EVERETT  MA  02149-5900 


1  BATTELLE 

TWSTIAC 
HUGGINS 
505  KING  AVE 
COLUMBUS,  OH  43201-2693 

1  COHEN  PROFESSIONAL  SERVICES 

ATTN  N  S  COHEN 
141  CHANNING  ST 
REDLANDS  CA  92373 

1  EXXON  RSCH  &  ENGRG  CO 

ATTN  A  DEAN 
RTE22E 

ANN  AND  ALE  NJ  08801 

1  GENERAL  APPLIED  SCIENCE  LAB  INC 

77  RAYNOR  AVE 
RONKONKAMA  NY  11779-6649 

1  GENERAL  ELECTRIC  ORD  SYSTEMS 

ATTN  JMANDZY 
100  PLASTICS  AVE 
PITTSFIELD  MA  01203 

1  GENERAL  MOTORS  RSCH  LAB 
ATTN  T  SLOANE 
PHYSICAL  CHEMISTRY  DEPT 
WARREN  MI  48090-9055 

2  HERCULES  INC 
ATTN  W  B  WALKUP 
E A  YOUNT 

ALLEGHENY  BALLISTICS  LAB 
ROCKET  CTR  WV  26726 

1  HERCULES  INC 

ATTN  R  V  CARTWRIGHT 
100  HOWARD  BLVD 
KENVIL  NJ  07847 

1  ALLIANT  TECHS YSTEMS  INC 

ATTN  D  E  BRODEN 
MS  MN50  2000 
MARINE  SYSTEMS  GROUP 
600  2ND  ST  NE 
HOPKINS  MN  55343 
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ALLIANT  TECHS  YSTEMS  INC 
ATTN  R  E  TOMPKINS 
MN  11  2720 
600  SECOND  ST  N 
HOPKINS  MN  55343 

IBM  CORPORATION 
ATTN  A  C  TAM 
RSRCH  DIV 
5600  COTTLE  RD 
SAN  JOSE  CA  95193 

IIT  RSRCH  INSTITUTE 
ATTN  R  F  REMALY 
10  W  35TH  ST 
CHICAGO  E  60616 

LOCKHEED  MISSEES  &  SPACE  CO 
ATTN  GEORGE  LO 
DEPT  52  35  B204  2 
3251  HANOVER  ST 
PALO  ALTO  CA  94304 

OLIN  ORDNANCE 
ATTN  V  MCDONALD  LIB 
PO  BOX  222 

ST  MARKS  FL  32355-0222 

PAUL  GOUGH  ASSOCIATES  INC 
ATTN  P  S  GOUGH 
1048  S  ST 

PORTSMOUTH  NH  03801-5423 

HUGHES  AIRCRAFT  CO 
ATTN  T  E  WARD 
8433  FALLBROOK  AVE 
CANOGA  PK  CA  91303 

ROCKWELL  INTERNATIONAL  CORP 
ATTN  J  E  FLANAGAN  HB02 
ROCKETDYNE  DIV 
6633  CANOGA  AVE 
CANOGA  PK  CA  91304 

SCIENCE  APPLICATIONS  INC 
ATTN  R  B  EDELMAN 
23146  CUMORAH  CREST 
WOODLAND  HELS  CA  91364 


3  SRI  INTERNATIONAL 
ATTN  G  SMITH 
D  CROSLEY 
D GOLDEN 

333  RAVENSWOOD  AVE 
MENLO  PK  CA  94025 

1  STEVENS  INSTITUTE  OF  TECH 

ATTN  R  MCALEVY  m 
DAVIDSON  LAB 
HOBOKEN  NJ  07030 

1  SVERDRUP  TECH  INC 

ATTN  R  J  LOCKE  MS  SVR  2 
LERC  GROUP 
2001  AEROSPACE  PKWY 
BROOK  PK  OH  44142 

1  SVERDRUP  TECH  INC 

ATTN  J  DEUR 
2001  AEROSPACE  PKWY 
BROOK  PK  OH  44142 

3  THIOKOL  CORPORATION 

ATTN  R  BIDDLE 
R  WELER 
TECH  LIB 
ELKTON  DIV 
PO  BOX  241 
ELKTON  MD  21921 

3  THIOKOL  CORP 

ATTN  S  J  BENNETT 
WASATCH  DIV 
PO  BOX  524 

BRIGHAM  CITY  UT  84302 

1  UNITED  TECH  RSCH  CTR 

ATTN  A  C  ECKBRETH 
E  HARTFORD  CT  06108 

1  UNITED  TECH  CORP 

ATTN  R  R  MILLER 
CHEMICAL  SYSTEM  DIV 
PO  BOX  49028 
SAN  JOSE  CA  95161-9028 

1  UNIVERSAL  PROPULSION  CO 

ATTN  H  J  MCSPADDEN 
25401  N  CENTRAL  AVE 
PHOENIX  AZ  85027-7837 
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1  VERITAY  TECH  INC 

ATTN  E  B  FISHER 
4845  MILLERSPORT  HWY 
E  AMHERST  NY  14051-0305 

1  FREEDMAN  ASSOCIATES 

ATTN  E  FREEDMAN 
2411  DIANA  RD 
BALTIMORE  MD  21209-1525 

6  ALLIANT  TECHS YSTEMS 

ATTN  J  BODE 
C  CANDLAND 
L  OSGOOD 
R  BURETTA 
R  BECKER 
M  SWENSON 
600  SECOND  ST  NE 
HOPKINS  MN  55343 

1  US  ARMY  BENET  LAB 

ATTN  SAM  SOPOK 
SMCAR  CCB  B 
WATERVLIET  NY  12189 


ABERDEEN  PROVING  GROUND 


36  DIR  USARL 

ATTN  AMSRL  WTP.AW  HORST 
AMSRL-WT-PC, 

R  A  FIFER 

GF  ADAMS 

W  R  ANDERSON 

R  A  BEYER 

S  WBUNTE 

C  F  CHABALOWSKI 

K  P  MCNEILL-BOONSTOPPEL 

A  COHEN 

R  CUMPTON 

R  DANIEL 

D  DEVYNCK 

N  F  FELL 

B  E  FORCH 

J  M  HEIMERL 

A  J  KOTLAR 

MRMANAA 

W  F  MCBRATNEY 

K  L  MCNESBY 

S  VMEDLIN 

M  S  MILLER 

A  W  MIZIOLEK 

S  H  MODIANO 

J  B  MORRIS 

J  E  NEWBERRY 

S  A  NEWTON 

R PATEL 

R  A  PESCE  RODRIGUEZ 
B  M  RICE 
R  C  SAUSA 
M  A  SCHROEDER 
J  A  VANDERHOFF 
M  WENS1NG 
A  WHREN 
J  M  WIDDER 
C  WILLIAMSON 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ART,  Report  Number  AKL-TR-695 _ Date  of  Report  February  1995 _ 

2.  Date  Report  Received _ _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ _ _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.) _ _ _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate.  _ _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


